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We extend the geometric measure of quantum discord, introduced and computed for 
two-qubit states inE to quantify non-classical correlations in composite Gaussian states 
of continuous variable systems. We lay the formalism for the evaluation of a Gaussian 
geometric discord in two-mode Gaussian states, and present explicit formulas for the 
class of two-mode squeezed thermal states. In such a case, under physical constraints of 
bounded mean energy, geometric discord is shown to admit upper and lower bounds for a 
fixed value of the conventional (entropic) quantum discord. We finally discuss alternative 
geometric approaches to quantify Gaussian quadrature correlations. 
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1. Introduction 

Gaussian states of quantised fields and matterlike systems constitute a tiny class 
in the arena of general continuous variable (CV) quantum states, yet their role is 
recognised as a leading one in quantum technology applications 2 . They arc defined 
as those states with a positive-anywhere, Gaussian-shaped Wigner distribution in 
phase space. As positivity of the Wigner function is sometimes regarded as a sig- 
nature of "classicality"^, one might wonder how come quantum protocols can be 
implemented and enhanced by the use of states which do not appear to radically 
deviate from a classical description. However, this picture is too simplistic. A more 
refined notion of classicality can be adopted, whereas quantum states p are con- 
sidered as "classical" if their Glauber-Sudarshan P representation, implicitly de- 
fined by p = J d 2 aP(a)\a}(a\ (with \a) denoting a coherent state), is regular and 
nonnegative^. With respect to this criterion, a numerical study has demonstrated 
that essentially almost all two-mode Gaussian states are non-classicaP. The typical 
example of a non-classical Gaussian state is a squeezed state. Squeezing represents 
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an essential ingredient to establish entanglement between two or more modes in a 
global Gaussian stated entanglement, in turn, constitutes an unambiguous quan- 
tum trait, ultimately enabling better-than-classical realisations of communication 
primitives such as quantum teleportatiorPl For this reason, the structure and dis- 
tribution of bipartite and multipartite entanglement in multimode Gaussian states, 
and its operational significance for teleportation and related schemes, has been 
extensively studied in recent years^. 

One could then ask whether non-entangled Gaussian states are just "classical" . 
Surprisingly enough, this question was settled only very recently. The concurrent 
works of Giorda and PariiPand Adesso and DattcP^ provide a clear, negative answer 
to that. In those studies, a paradigm more general than the entanglement-versus- 
separability one^2 has been explored to assess the quantumness of correlations in 
bipartite Gaussian states. This paradigm was proposed one decade ago in two pio- 
neering independent works by Zurek and Vedral (and their collaborators pHEH They 
identified signatures of quantumness stemming directly from the non-commutativity 
of quantum observables, which can be retrieved for instance in the fact that local 
measurements on one or more subsystems of a quantum composite system typically 
induce a disturbance on the state of the system, at striking variance with what 
happens when extracting marginal information from a classical joint probability 
distribution. This generally happens in entangled as well as in generic separable 
quantum states. Precisely, this is the case for almost all quantum states, and the 
purely classically correlated (or "classical") states constitute just a null measure, 
nowhere dense subset of any composite Hilbert spaced The classical states are just 
those with a 'collapsed' wavefunction, also known as pointer states^, whose density 
matrix is diagonal in some local product basis. For a more precise mathematical 
definition, please refer to the paper by Gharibian et al. in this special issued 

The evidence of quantum correlations in separable states has been mostly ig- 
nored until it was speculated that such correlations might be at the root of the 
quantum speed-up in computational algorithms that run over mixed qubit regis- 
ters, such as the DQC1 modeP^. In an attempt to capture, operationally and/or 
from a more abstract perspective, the essence of general non-classical correlations, 
several measures have been proposed for their quantification, among which we men- 
tion Refs. I 12 ! 13 ! 18 ! 19 ! 20 ! 21 ! 22 ! 1 ! 2 ^ F r two-mode Gaussian states, detailed results are 
available concerning the evaluation of the quantum discord and of the (ame- 
liorated) measurement-induced disturbance^. In the case of quantum discord T> 
E^l, which is one of the most popular measures of non-classical correlations at the 
moment, a closed formula is available for general Gaussian states of two-mode CV 
systems, under the restriction of Gaussian local measurements performed on one 
of the partjff^l, while only a semi-analytical expression exists for the quantum dis- 
cord of general states of two qubits^. On the other hand, an intuitive geometric 
approach to the quantification of non-classical correlations has recently led to the 
introduction of the "geometric measure of quantum discord" (or geometric discord) 
Dq I 1 ! 26 !, which is easily computable in closed form for general two-qubit states. 
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In this paper we propose a generalisation of the geometric discord to Gaussian 
states of CV systems. We define the "Gaussian geometric discord" (GGD) as the 
minimum squared Hilbert-Schmidt distance between a Gaussian state and the clos- 
est 'classical-quantum' state obtained after a local generalised Gaussian positive- 
operator-valued-measurement (GPOVM) performed on one party only. Focusing on 
two-mode Gaussian states, we provide a covariance matrix formulation of the prob- 
lem, and we solve the optimisation explicitly for the relevant subclass of two-mode 
squeezed thermal states: the optimal measurement in this case is shown to be a het- 
erodyne detection, similarly to the case of conventional quantum discordlES, and 
the GGD turns out to be symmetric under swapping of the two parties. We then 
compare GGD with the quantum discord for two-mode squeezed thermal states. 
While an upper bound on GGD can be identified exactly at fixed quantum discord, 
states can be found with arbitrarily small GGD for any arbitrarily large quantum 
discord. We recognise this as a consequence of the known result that infinitely 
entangled states of CV systems are dense in norm in infinite-dimensional Hilbcrt 
spaces^. In fact, imposing simple physical constraints of finiteness of the mean 
energy, we are able to recover as well a hierarchy of lower bounds on GGD ver- 
sus discord. We finally discuss possible alternative generalisations of the geometric 
discorcGEH to Gaussian states. If one applies, for instance, the notion of minimum 
Hilbert-Schmidt distance from classical states directly at the level of covariance ma- 
trices, an easily computable measure of total quadrature correlations is obtained. 

The paper is structured as follows. In section [2] we recall the notion of geomet- 
ric discord, and propose its generalisation to Gaussian states of CV systems. In 
section [3] we evaluate the GGD for two-mode Gaussian states, and compare it to 
the conventional (entropic) quantum discord. In section [4] we discuss an alterna- 
tive geometric measure of Gaussian correlations. Finally, section [5] contains a brief 
summary and outlook. 

2. Geometric discord 

The geometric measure of quantum discord Dq has been recently proposed as a 
simple and intuitive quantifier of general non-classical correlations in bipartite quan- 
tum states^. Let us consider a bipartite system AB in the state pab and imagine 
to perform a local measurement Hb on B. Almost all (entangled or separable) 
states pab will be disturbed by any such measurement^. However, for the class of 
so-called classical-quantum state£lE51 of the form 



where pi is a probability distribution and {\i)} is a basis for the Hilbert space 
of subsystem B, there exists at least one local measurement that leaves the state 
invariant. Denoting by A be the set of classical-quantum states with respect to the 
A versus B bipartition, the geometric discord Dq is defined as the squared Hilbert- 
Schmidt distance between the state pab and the closest classical-quantum state 




(1) 
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XAB®, 



Dg(pab) = inf \\pab - Xab\\ 2 ■ ( 2 ) 

XAB&A 



where ||M|| 2 = ^tv(MM^) stands for the 2-norm, or Hilbert-Schmidt norm, of 
the matrix M. The quantity Dq(pab) in Eq. ^ vanishes on classical-quantum 
states and is upper bounded by the state purity tx(p\ B ) in general. Interestingly, 
the geometric discord can be exactly reformulated as the minimal disturbance, 
measured in terms of the squared Hilbert-Schmidt distance, induced on the state 
Pab by any projective measurement II s on subsystem B S3 

Dg (pab )= inf || pab - IT B (pab) || \ ■ (3) 

We notice that the geometric discord is in general not symmetric under the swap 
of the two parties, A -h- B, i.e. measuring A rather than B may induce different 
amounts of disturbance on generic bipartite states. 

In this paper we shall consider pab to be a bipartite, ro-mode Gaussian state 
Pab> i- e - a state with Gaussian Wigner distribution on the quantum phase space 
of n quantised harmonic oscillators, partitioned into two blocks of ua and ns 
modes, respectively (with n = ua + Hb). Without loss of generality, we can assume 
that all the hrst moments of the canonical phase-space operators are set to zero, 
because we are aiming at computing correlations between the modes, which are not 
affected by local displacements. The states under consideration are then entirely 
described by their covariance matrix Sab of elements = ti(pAB{Ri, Rj}+) 
where R = (xi,pi, . . . , x n ,p n ) is the vector of phase-space operators satisfying the 
canonical commutation relations [Ri,Rj] = iClij, with Q = fjl)* being the 
symplectic matrix^. 

We then define the Gaussian geometric discord (GGD) of a Gaussian state 
Pab as m Eq- with TIb constrained to be a generalised GPOVM 11^ , i.e. a map 
sending Gaussian states into Gaussian states. These measurements coincide with the 
standard toolbox of linear optics, i.e., can be realised using beam splitters, phase 
shifters, squeezers, appending ancillary vacuum modes, and performing balanced 
homodyne detection^!. The GGD is then, explicitly, 

D G (p%)=ird\\p s AB -IL%(p AB )f 2 , (4) 

which is normalised between and 1 for Gaussian CV states. From now on, we 
drop all the superscripts "Q" and always assume, implicitly, that all the states and 
measurements we consider are Gaussian. 

3. Gaussian geometric discord for two-mode Gaussian states 

We shall focus on two-mode Gaussian states pab , specified by their 4x4 covariance 
matrix 
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which has to be positive definite and has to comply with the uncertainty prin- 
ciple, Sab + ifi > 0, to ensure its correspondence to a physical state. The 
GPOVM on the single- mode subsystem B can be written in general as Hb(v) = 
tt~ 1 Wb{j])wbW\{vi) where Wb(i]) — exp(ry6^ — r)*b) is the Weyl operator, b = 
(£b + ipB)/V%, 7T _1 / rf 2 ijng(i/) = 1, and the seed tub is the density matrix of a 
(generally mixed) single-mode Gaussian state with zero mean and covariance ma- 
trix <j b- After the measurement, the state of mode B is projected onto the seed wb, 
while the posterior state w a of mode A has a covariance matrix that is crucially not 
dependent on the measurement outcome^, being given by the Schur complement 

<?A = ol - j(/3 + cr B y l l T • (6) 

This entails that no classically correlated mixture can be realised as the outcome of 
a local Gaussian measurement, and the only classical-quantum two- mode Gaussian 
states are just uncorrelated tensor product states, as already proven irP^l. The GGD 
for two-mode Gaussian states thus reads 

Dg{pab) = inf \\pab — ^A ® ^b\\\ , (7) 

where the single-mode states wa,b are introduced above. 

We remark that in our notation, the optimisation is over tub only, and there 
is no free parameter left in wa- Let us comment on this issue in some detail. 
One might think to define alternatively the Gaussian counterpart of the geometric 
discord as in Eq. ([2]), by letting xab be any possible two- mode tensor product 
Gaussian state, i.e. xab = Xa ® Xb, with completely arbitrary xa and Xb, and 
optimising over both. This class of states is clearly more general than the subset 
of tensor product states emerging from the GPOVM analysis, which enters in our 
definition Eq. ([7]), thus possibly enabling a further minimisation of the Hilbert- 
Schmidt distance compared with the optimal value of Dq obtained from Eq. 
However, a numerical comparison [see Fig. [I] reveals that the two definitions re- 
sult in actually quite close (although not exactly coincident) values for randomly 
generated two-mode Gaussian states. Therefore, we believe it more natural and 
operationally wise to stick with the definition of Dq that we introduced in Eq. Q. 
A reconcilation between Eqs. ^ and ([3| for bipartite Gaussian states would be 
possible only enlarging the set of operations Hb to general, even non-Gaussian CV 
measurements^ (and, equivalently, the set A to generally non-Gaussian two-mode 
classical-quantum states), which is beyond the scope of this paper. 

We now proceed with the evaluation of Eq. ([7]). Recalling the formula for the 
overlap of two Gaussian states pi and p2 with covariance matrices cr\ and <r 2 ^ 
ti(p 1 p 2 ) = l/ v /det[(<r 1 + <r 2 )/2], one has 



D g {-£ab) = inf • 



[ Vdet Tsab y/det(a A ©<x s ) ^/det[(S AB + cr A © <tb)/2] 
with Hab given by Eq. Q and cta given by Eq. Q. 



(8) 
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Fig. 1. Comparison between two possible Gaussian extensions of the geometric discord (see the 
main text for further details), one emerging from Eq. |2| (vertical axis), and the other emerging 
from Eq. (J3J (horizontal axis). The latter is what we define as GDD in Eq. fil and adopt in this 
paper. The points correspond to 10000 randomly generated two-mode Gaussian states. All the 
quantities plotted are dimensionless. 



The measure Dq is local-unitarily-invariant by definition. On the other hand, by 
means of local unitary operations (symplectic transformations in phase space) , the 
covariance matrix of every two-mode Gaussian state can be brought into a unique 
standard forrrP^, 



/a0c0\ 
OaOd 
c b 

\0d0bj 



(9) 



with a, b > 1, y/ab-1 >c> \d\. There is no loss of generality in restricting ourselves 
to states in standard form in order to proceed with the explicit computation of 
Dq. The formulas that we derive in the following are in fact valid for arbitrary 
two-mode Gaussian states, provided one identifies the standard form covariances 
in terms of the four symplectic invariants det a = a 2 , det/3 = b 2 , det7 = cd, 
det S AB = (ab - c 2 )(ab - d 2 )®. 

Now, the seed state is in general a single-mode rotated, squeezed, thermal state 
with covariance matrix 

( m XcOS 2 (9) + ^^ _ m(A'-l)co BW Bin W 

aB \ _"(* 3 -iW«) ■*-(«> r^ +TOA sin 2 (tf) 1 ' 

with m > 1 and A > 0. The task we are facing is then the minimisation of Eq. ([8j for 
two-mode states in standard form, over the GPOVM parameters m (temperature), 
A (squeezing), and 9 (rotation angle). The optimisation over 9 is straightforward 
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and yields 9 — 0. Eq. (18]) then becomes 



D G {H AB ) = inf 

m,A 



m 2 (a(b+\m)-c 2 )(a(b\+m)-d 2 \) I (a(b+\m) — c 2 )(a(b\+m) — d 2 \) 

(b\+m)(b+\m) \ X 

(11) 



y/(ab-c 2 )(ab-d 2 ) 

The remaining optimisations can be in principle solved analytically, although 
the resulting formulas are not very handy in general and there is no need to report 
them here explicitely. We prefer instead to specify our attention to the physically 
relevant subclass of two-mode squeezed thermal states, characterised by d = ±c 
in their covariance matrix For these states, the problem admits a simple 

solution, and the GGD turns out to be minimised at 



Vab(y/4ab - 3c 2 + sfab) 

A = 1 , m = . 

3a 

The least disturbing GPOVM for squeezed thermal states, according to the Hilbert- 
Schmidt distance, is thus a (noisy) heterodyne detection, a result which is analogous 
to what found in the case of quantum discordH112l The GGD of two-mode squeezed 
thermal states finally acquires the following compact expression 

= ^bh - ( , 9 r-V • (12) 

U/Aab - 3c 2 + Vabj 

Notice that the case of pure two-mode squeezed states is recovered for b = a,c = 



\Jo? — 1. A quite unexpected (yet nice) feature of Eq. ( 12 1 is that Dq of two-mode 
squeezed thermal states happens to be symmetric under party swap [a O b), despite 
the generally nonsymmetric definition of the GGD. 

We can test the accuracy of the GGD in quantifying non-classical correlations of 
Gaussian states by means of a comparison with the conventional quantum discord 
T>^, whose definition and evaluation for two-mode Gaussian states can be found 
in detail iri 9 l 10 [ Here, it suffices to recall that the (Gaussian) quantum discord T> 
quantifies the minimal change in the mutual information of a bipartite (Gaussian) 
state after an optimal, least disturbing local (Gaussian) measurement is performed 
on subsystem B; notice that T> can range up to infinity, while Dq is normalised to 
one as already remarked. We run the analysis for squeezed thermal states, although 
the conclusions can be easily extended to more general two-mode Gaussian states. 

A first observation that one can make, is that there is no a priori finite lower 
bound on Dq for a fixed amount of quantum discord T>. This can be understood by 
recalling that in infinite-dimensional Hilbert spaces, infinitely entangled states (or, 
equivalently, states with diverging quantum discord^) constitute a (trace-norm) 
dense set. Our results show that the same holds for the Hilbert-Schmidt norm as 
well. This means, in particular, that in any neighbourhood of every product (clas- 
sical) state xab lies an arbitrarily strongly entangled (non-classical) state pab- 
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The latter state pab would thus have arbitrarily large quantum discord, and ar- 
bitrarily small geometric discord. We can provide an example family of two-mode 
squeezed thermal states that falls into this particular category. Consider the states 
with standard form covariances 



b=l + e, c= y/(a+l)e, with < e < a - 1 . (13) 



A plot of Dq versus V is shown in Fig. 2(a) for the states of Eq. (13 1, with a = 2 



and k = 1,2,..., 10 (from left to right, solid lines). For these states, T>(a,e) 



4(a+l)tanh 1 ( 2a _ £ e+2 ) +2e(-a+e+l) coth 1 (a-e)-e(e+2) log (^j 



and Da(a,e) can be obtained from Eq. (12). By setting e = x(a — 1) with a finite 



/[2(2+e)] 



ratio < x < 1, one sees that, rigorously Dq — > for this class of Gaussian 

a— >oo 

states, and T> — > —x^ 1 ln(l — x), which increases arbitrarily for x approaching 1, 

a— >oo 

thus demonstrating the claim. Let us remark that this feature is general, only due 
to the geometry of CV state spaces, and is not a peculiarity of our definition Eq. Q 
of Do- Even allowing for any non-Gaussian measurement on B, or equivalently any 
non-Gaussian classical-quantum state in Eq. ([2| , one could only further reduce the 
geometric discord, which would then still vanish for all those cases, like the two- 
mode family considered above, in which the GGD defined by Eq. ^ vanishes in 
the presence of an arbitrary degree of non-classical correlations measured by T> (or 
by some entanglement measure). 

A physically motivated way to address this issue is by introducing a cap on the 
mean energy per mode of the Gaussian states under investigational. This corre- 
sponds to imposing tr(<r) < 2N on the marginal covariance matrix of any single- 
mode Gaussian state, for some finite N > 1. The pathological behaviour in fact oc- 
curs in the unphysical limit of diverging mean energy. For squeezed thermal states 
in standard form, the finite energy constraint simply imposes a,b < N. We have 
run an extensive numerical comparison of Dq versus T> for 10 6 randomly generated 
two-mode squeezed thermal states with bounded mean energy per mode (setting 



N = 25), as shown in Fig. 2(b) In this expected, a well definite region in 

the parameter space is filled by the random states, and an exact lower bound on 
the GGD can be identified at fixed quantum discord. Interestingly, it is saturated 



precisely by states of the form Eq. ( 13 ) with a = N. 

We have also observed the presence of a general upper bound on Dq versus T> 
for two-mode squeezed thermal states. It does not depend on the energy capping 
(i.e., it is independent of N) and thus stands as an upper limit on the GGD for any 
two-mode squeezed thermal state, even asymptotically infinitely entangled. The 
states saturating the upper bound can be found within the two-parameter class 
characterised by Eq. (|9| with 

a = b-e, c= y/(b+l)(b-e-\), with < e < b - 1 . (14) 

For these states, V(b,e) = §[(1 - b) ln(6 - 1) + {b + l)ln(6 + 1) + eln(e) - (e + 



2) ln(e + 2)] and Dc{b,e) can be found again via Eq. (12). We are interested in 
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D 

(a) 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 

D 

(b) 



Fig. 2. Upper (dotted black) and lower (solid black) bounds on the GDD Dq at fixed quantum 
discord T> for two-mode squeezed thermal Gaussian states with (a) mean energy per mode bounded 

by N = 2 fc with k = 1, 2, , 10 (solid lines, from left to right) and (b) mean energy per mode 

bounded by N = 25. In both plots, the lower boundaries (solid lines) correspond to states of the 
form Eq. \13\ , while the upper boundary (dotted, independent of N) correspond to optimised 
states belonging to the class of Eq. | |14| , as detailed in the main text. Panel (b) includes a 
numerical exploration of Dq versus T> for 10 6 randomly generated two-mode squeezed thermal 
states with mean energy per mode bounded by TV = 25 (gray points). All the quantities plotted 
are dimensionless. 

a one-parameter subclass of this family, for which Dq is maximised at a fixed 
value of T>. This optimisation can be performed numerically in two steps: first, 
identifying the contours of constant discord T> in the space of parameters b and e; 
then, maximising Dq along each contour. The result is plotted as a dotted black 
curve in Fig. [2} We notice that, for P>0, the maximal Dq saturates quite rapidly 
to the maximal value of 1, and the states sitting on the upper boundary converge 
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simply to pure two-mode squeezed states, i.e., e = in Eq. (14 1 is optimal for 
highly correlated states. Again, this general result would still hold if the definition 
m Eq. (2| was employed. 

Summarising, if no bound is imposed on the mean energy per mode, what 
we have found is that, essentially, for strongly non-classically correlated Gaussian 
states, the GGD is a rather insensitive measure of quantum correlations, as it can 
range arbitrarily between and 1. On the other hand, for finite degrees of cor- 
relations (on entangled as well as separable states), respecting a physical cap on 
the mean energy, nontrivial upper and lower bounds exist on the GGD at fixed 
quantum discord, as demonstrated explicitly for the family of two-mode squeezed 
thermal states. We remark that exact upper and lower bounds on the geometric 
discord at fixed quantum discord have been recently demonstrated in the case of 
arbitrary states of two-qubit systems^. 

4. Geometric quadrature correlations 

In this section we briefly discuss an alternative possible pathway to extend the 
notion of a geometric measure of non-classical correlations to bipartite Gaussian 
states. One might apply the notion of minimal Hilbert-Schmidt distance from the set 
of classical-quantum states, as in Eq. ([2]), directly at the level of covariance matrices, 
as they provide a compact and complete characterisation of general Gaussian states 
up to local unitaries. We then define the "geometric quadrature correlations" (GQC) 
as 

Qg(£ab) =inf JV AB -T AB \\ 2 2 . (15) 

where T denotes the set of all physical covariance matrices Tab associated to 
classical-quantum Gaussian states with respect to the A versus B bipartition of 
modes. This set can be characterised following the results of Ref. ^ (see also the 
Supplemental Material irP^). The quantity Qc in Eq. (15) is always nonnegative 



and manifestly invariant under local symplectic transformations, it vanishes on 
classical-quantum Gaussian states, and can range up to infinity. 

For two-mode Gaussian states, Qg simply reduces to the minimum distance, 
in phase space, from completely uncorrelated tensor product states: QgC^ab) = 
inf CTl , CT2 ||Eab — {cr\ ® f 2 )||2 ■ As in the case of the GGD, we can evaluate this 
quantity explicitly assuming Sab to be in standard form, Eq. (|9|), and optimising 



the distance over arbitrary single-mode states <t 12 of the form Eq. (10), specified 



by the sets of parameters {mi, Ai, 9\} and {rri2, A 2 , #2}, respectively. We then have 
Q g (£ab) = inf tr [(Sab - a x ® er 2 ) 2 ] . (16) 

{mi, 2, Ai, 2,61,2} 

Explicitly, QgCZab) = inf {mi 2 , Al aA 2> [ - 26A 2 (\ 2 + l) A?m + (A| + l) \\n\ + 
X 2 (2Af (a 2 + b 2 + c 2 + d 2 ) - 2a\ 1 (X 2 + l) m a + (A? + l) mf) ] /(A? A|). 

The expression is much more tractable than Dg and there is no dependence 
on the angles 9j . The remaining optimisation over rrij , A j can be easily solved by 



October 13, 2011 0:18 WSPC/INSTRUCTION FILE gaussgeomq 



Gaussian Geometric Discord 11 

imposing vanishing partial derivatives with respect to the parameters, and one 
finally gets, for arbitrary two-mode Gaussian states in standard form, that Qq 
attains its global minimum for Ai = A2 = 1, mi — a, rri2 = b. In other words, the 
closest classical state is just the tensor product of the two marginals, Tab =ct®/3. 
The GQC thus acquires a very simple form: 

QgCEab) = 2(c 2 + d 2 ) (17) 
_ 2(det a + det 2 7 - dct T, AB ) 
y/det a dct (3 

where in the second line we have explicitly rewritten the solution in terms of the 
invariants of the covariance matrix [Eq. j5b], in order to provide a general formula 
not relying on the standard form. 

The GQC is then a symmetric measure that obviously quantifies total correla- 



tions of two-mode Gaussian states, as intuitively clear from the definition Eq. (16 1. 
Its standard form expression is particularly transparent, as it just coincides with the 
sum of the squared off-diagonal quadrature correlations in the state characterised 
by Sab, which essentially account for all the possible correlations of Gaussian 
states (whose higher order correlations in the canonical phase-space operators are 
entirely specified by the ones between second moments). We notice a similarity 
between this formula and its counterpart applied at the level of density matrices 
for finite-dimensional systems, where quantum correlations can be quantified by the 
so-called "minimum entanglement potential" (or negativity of quantumness) , which 
amounts to the sum of the modulus square of all the off-diagonal coherences in the 
density matrix of a bipartite system, minimised over all possible local basest 

We clarify with an example how GQC in the Gaussian setting quantifies the 
amount of total (i.e. both classical and quantum) correlations. Let us consider the 
family of two-mode Gaussian states with standard form covariances given by 

b = a, c = a-l, d = 0. (18) 

For these states, both the quantum discord^ T> and the GGD Dq [Eq. ([8])] arc 
very small for any a > 1, never exceeding 1/40 and vanishing exactly in the limit 
a — > 00. On the contrary, the GQC quantity Qq scales as 2(a — l) 2 , diverging for 
large a. The measure Qq is then closer in spirit to the "bit quadrature correlations" 
studied iiP^l for two-mode CV states. 



5. Conclusions 

In this paper we addressed the quantification of general non-classical correlations 
in Gaussian states of CV systems from a geometric perspective. Generalising the 
known finite-dimensional definitions' 1 ! 26 !, we proposed a Gaussian version of the 
geometric measure of discord, defined as the minimum distance between a bipartite 
Gaussian state and the closest classical-quantum Gaussian state obtained after a 
local Gaussian measurement on one subsystem only. We calculated the measure 
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for two-mode Gaussian states, focusing in particular on squeezed thermal states. 
We then compared the Gaussian geometric discord with the conventional Gaussian 
quantum discord-!- 2 ! 9 ! 10 !, and identified upper and lower bounds on one quantity as 
a function of the other, upon constraining the mean energy per mode. In absence 
of such a provision, the geometric discord can vanish for arbitrarily strongly non- 
classical (entangled) Gaussian states, as a consequence of the geometry of CV state 
spaces^. 

We finally analysed a different Gaussian counterpart to the geometric discord, 
where the distance from the set of classical-quantum states is calculated directly 
on covariance matrices; this operation results in a simple quantifier that captures 
the total quadrature correlations rather than the genuine non-classical correlations 
of Gaussian states. 

Contrarily to the two-qubit case, where the geometric quantum discord is easier 
to calculate than the conventional discordll 26 ! 25 !, we find no significant advantage 
(beyond the fact of managing polynomial rather than logarithmic functions) in 
choosing a geometric perspective rather than an entropic one, for the quantifica- 
tion of Gaussian non-classical correlations. It would be interesting to investigate al- 
ternative distance-based measures of non-classical correlations for Gaussian states, 
based on different norms, e.g. in terms of phase-space distinguishability measures. 

We finally mention that, following the analysis o f 1 ! 33 !, quantum discord has 
been very recently witnessed in qubit setups implemented with nuclear magnetic 
resonance^. Inspired by that, an interesting future task could be trying to devise 
proper quantitative witnesses of Gaussian (geometric) quantum correlations, ac- 
cessible experimentally in all-optical implementations or hybrid setups interfacing 
light modes and ultracold atomic ensembles. 
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